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Abstract 

Inhomogeneous and anisotropic cosmologies are modeled withing the frame- 
work of scalar-tensor gravity theories. The inhomogeneities are calculated 
to third-order in the so-called long-wavelength iteration scheme. We write 
the solutions for general scalar coupling and discuss what happens to the 
third-order terms when the scalar-tensor solution approaches at first-order 
the general relativistic one. We work out in some detail the case of Brans- 
Dicke coupling and determine the conditions for which the anisotropy and 
inhomogeneity decay as time increases. The matter is taken to be that of 
perfect fluid with a barotropic equation of state. 



PACS number(s): 04.50.+h, 98.80.H 
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I. INTRODUCTION 



Recently, Damour and Nordvedt |TJ have shown quite generally that General Relativity 
is an attractor of a large class of scalar-tensor gravity theories in the context of cosmology 
(see also ||). Earlier, Bekenstein and Meisels || had shown this to be true for Bekenstein's 
variable mass theory |4| . In both studies only Friedmann- Robertson- Walker spacetimes were 
considered. We shall here extend these results to anisotropic and inhomogeneous universes 
using the long-wavelength iteration scheme. A general barotropic equation of state is used 
so that universes driven by ordinary as well as "inflationary" matter can be considered. 

The long-wavelength iteration scheme is an approach that allows for a generic description 
of cosmological inhomogeneities and anisotropies that goes beyond the linear regime. This 
approach was introduced by Lifschitz and Khalatnikov H] and Tomita [0]. Recently various 
groups have employed what amounts to the same scheme: Salopek, Stewart and collaborators 
(see J/J and references therein) use it in their analysis of the Hamilton- Jacobi equation 
for general relativity; Comer, Deruelle, Langlois and collaborators (see J5],|9| and references 
therein, as well as [|lOj and discuss among other things how the equation of state for 



matter affects the growth or decay of cosmological inhomogeneities (see also [12-T^]); finally 



Soda et al. [pl| solve the Hamilton- Jacobi equation in Brans-Dicke theory when matter is 
pressureless. We shall here extend the Soda et al. result to fluids with pressure and to other 
couplings, using the approach developped in [|J and |§. 

The essential point of the long-wavelength iteration scheme is that spatial gradients must 
be small in the following sense: Choose a synchronous gauge where the metric takes the form 

ds 2 = -dt 2 + 7ij(t, x k )dx L dx j . (1) 

Let L be the characteristic comoving length on which the spatial metric 7^ varies spatially: 
di^ik ~ L~ l ^jk. The Hubble time H^ 1 (where H = a/a, a 2 = 7 1//3 , 7 being the determinant 
of 7^ and an overdot signifying d/dt) is the characteristic time on which the metric varies 
temporally: 7^ ~ Hjij. The long- wavelength approximation, accurate when L ^> Rh 
where Rh = (aif ) _1 is the local (comoving) Hubble radius, consists in neglecting all spatial 
gradients in the Einstein equations. However, the Hubble radius, in general, depends on 
time and eventually the condition L 3> Rh can be violated. To get more accurate (and 
"longer- lived" ) results, 7^ is expanded in terms of ever increasing spatial gradients of an 
arbitrary "seed" metric. Since each term in the expansion must be a tensor, the gradients 
of the "seed" enter via all possible contractions and derivatives (leaving two free indices) 
of the Riemann (or equivalently in 3-dimensions, Ricci) tensor formed from the "seed". 
The coefficients of the gradient terms entering the expansion are time- dependent and are 
determined order by order. For more discussion, see 

The only difficult part about implementing the iteration scheme is finding exact solutions 
at the very first iteration. Fortunately Barrow and Mimoso [16|] have given an algorithmic 
way of generating such solutions for scalar-tensor gravity theories. We shall adopt their 
techniques and extend them to the third-order. To get an explicit solution however the 
scalar-tensor coupling must be specified. We shall work out in some detail the Brans-Dicke 
coupling, consider the class studied by Damour and Nordvedt |l| and comment on those 
chosen by Barrow and Mimoso [|n| and Serna and Alimi [I7|. 

In Sec. we give the action for scalar-tensor gravity in the so-called "Jordan frame" and 
write the field equations that follow in a synchronous gauge. In Sec. [Hj we show how to solve 
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them in general in the long-wavelength approximation following Barrow and Mimoso [16 



and particularize to the specific couplings mentioned above. In Sec. [IV] the second iteration 
equations are written down and their solutions explored in the late time limit. Finally the 
corresponding calculations in the "Einstein frame" are outlined in an Appendix. 

The conclusion will be that, as in general relativity, inhomogeneities decay only if matter 
is inflationary. However there is an additional condition on the scalar-tensor coupling which 



is it must be close enough to its general relativistic value (see eq. in the text 



II. THE FIELD EQUATIONS 

In the "Jordan frame" the action for scalar-tensor gravity is 

S = ^r I d*xy/=g[<f>R - 0-^(0)^0^0] + S M ■ (2) 

1D7T J 

The function cu(0) is an a priori arbitrary function which couples the scalar field to the 
metric g^ v . Sm is the action for the matter which will be taken to be a perfect fluid. Its 
stress-energy tensor is 

T^ u = (p + p)u^u u + pg^ (3) 

where is its 4-velocity (u^u^ = —1) and where the energy density p and pressure p are 
chosen to be related by a barotropic equation of state of the form 

p = (r - i) P (4) 

where T is a constant (0 < T < 2). 

The field equations which follow from extremising eq. (0) are 

R " = j ( T » - \ T ^) + \ {€ + + ^ W (5) 

and 

fi 1/2 (ft 1/2 cT0) ;(T = 8ttT (6) 

where Q = 2u + 3 (Q is assumed to be positive) and T is the trace of the stress-energy 
tensor. 

In the synchronous gauge (see eq. ([!])), letting 7 U (which raises spatial indices) be the 
inverse of 7^ (which lowers spatial indices), defining = 7^, K = K\, the field equations 
become 

\k + \k{K< = f (l - - Tu\) P - i (30 + \K* -4)- ^ , (7) 



X -{K\ - K8i) }j = jT^l + uiUjpUt + ^{Kl^j - 20|.) - ^00,* , (8) 
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-m- 



+ 2 K * - 4 



+ 



si 



(9) 



" Ft {Ql/2 ^ ~ \ K ^* + (fil/20 ' 4)|j = J^ (3r " 4)p (10) 

where ( 3 'R%j is the 3-Ricci tensor formed out of 7^, a slash means a covariant derivative with 
respect to 7^ and 7 = det(Ty). Eqs. (|7|-|T0D are 11 equations for the 11 unknowns p, Ui, 
and jij. 



III. THE LONG- WAVELENGTH FIELD EQUATIONS AND THEIR SOLUTIONS 

In the long-wavelength approximation all spatial gradients in eqs. (|7|-|i~0D are ignored. 
Eq. (^|) then simply implies that 



m = . 



The trace-free part of eq. (0) reduces to 



dt 



Vj<j>[Kl- -K5j 



(11) 



(12) 



which yields 



A A 1 

—5- H 

A A s / 2 c 



■S:> 



(13) 



where we have defined A = a 2 = 7 1//3 so that K = 3 A/ A, and where the so-called anisotropy 
matrix Sj is time-independent and traceless. Eq. ( |13"D can be integrated to yield the 3-metric 
jij (see for details). The final result is 



(14) 



where the ef, a = 1,2,3, are the components of the 3 eigenvectors of Sj and 7 a b = 
A Diag [Aie rir , A 2 e r2T , A 3 e r3T ] with r a the eigenvalues of Si (with r\ + r 2 + r 3 = because 
Si is traceless), A a three functions of space such that AiA 2 A 3 = 1, and dr/dt = (A 3 / 2 0) _1 . 
When S? = the result simplifies into 



A L 



(15) 



where hij(x k ) is an arbitrary "seed" metric. 

The solution for p is most easily obtained from the conservation law u u T£ = which 
follows from the field equations (|5|-|6]) and boils down, at lowest order in the gradients, to 
p + 3THp = (with H = A/2 A) and thus 
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(16) 



where Mi is a function of space only. Note that eq. flip]) holds whether SI is zero or not. 

What remains is to solve for A and 4>. They are given by the traces of eq. (|9]), eq. (^) 
and eq. fllPp, neglecting all gradients and using eqs. flT3] ) and fllBp. (One of the equations of 
course is redundant and can be used as a consistency check.) When Sj = these equations 
boil down to the Friedmann equations. Barrow and Mimoso |16| have shown how to reduce 
the problem to a quadrature when Q(t) rather than Q(<p) is specified. For Sj 7^ Barrow 
and Mimoso's procedure can still be followed the only difference being that their constraint 
equation acquires an additional anisotropy term. We therefore refer the reader to || for 
details. The final results are 

ln0=r^^- (17) 
Jo g{r}) 



and 



A J2 i n =M 1 (4-3T) 9 -^ (18) 





where rj is related to the time t by 



, gO-T) 



dr}/dt = A^^yJ3/n , (19) 

and where g(rj) is the function which characterizes the coupling through the relation 

n = 4-3r [2g> - (3r - 4)r/] 2 

3(2 -T) 2 4(7- (3r-4)r ? 2 + 5/S'j/[6M 1 2 (4-3r)] 1 J 

with g>' = dg/drj. Integration of this equation for g yields a constant of integration, which 
can be used to parametrize the space of solutions. The anisotropy matrix only shows up in 
eq. (|20|) . However, as one can see from this equation, the anisotropy will be negligible at 
late times. 



A. Brans-Dicke coupling 

This case corresponds to Q = const. The Friedmann solutions (Sj = in eq. (120|) ) have 
been extensively studied (see, e.g., |3|,[T8|-pll and references therein). From eq. (|20|) and the 



constancy of f2 we obtain the "generating" function g: 

4(4 - 3^(17) = [30(2 - T) 2 - (4 - 3r) V + *V + ^ _ r)2 " ^ (21) 

where A is a function of space which parameterizes the solutions. It is then an elementary 
exercise to integrate eqs. (|PTD and (p~8|) . Only the relationship between rj and t [eq. (|19|)1 
remains non trivial as it involves a non elementary function of the type / dr)(r)—r)+) a (r)—r)-)P 
where rj± are the roots of g(rj) = 0. At late time (or at all times if A = Sj — 0) the solution 
takes the familiar form: 
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x1 * q = 4 ( 4 - 3F ) ( 22 ) 

' q (2-3r)(4-3r) + 3r(2-r)fi ' 1 ; 



a oc h- 2[n(2-r)-(4- 3 r)] 
°" ' (2-3r)(4-3r) + 3r(2-r)fi • [ 6) 

(When T = 1, eqs. (|22]-|2"3|) agree with the solution discussed by Weinberg PIfl .) As A must 
remain positive we must have that (4 — 3T)g(rj) > (see eq. (|18|) ) and therefore 

Q > (4-3r) 2 /3(2-r) 2 . (24) 



(See [[22| for more discussion on the vanishing of anisotropy in Brans-Dicke theory.) 



B. Dynamical couplings 



Various coupling functions fi(0) have been considered in the literature (see e.g. [17 and 
references therein). To be viable the space of the corresponding cosmological solutions must 
exhibit a basin of attraction towards general relativity. This means that when, say, <fi — > 1 
then Q — > oo (with (dQ / d(f>) / fl 3 —> [p3] ). Many different behaviours of course are possible, 
for example 

Q^a/lncp (25) 
which is the one considered by Damour and Nordvedt ffl] and also by Barrow and Mimoso 



T6| as it is the late-time behaviour of all solutions generated by a function g oc r\ n with 



n > 2. As follows from eqs. (17-19) we have indeed [16 



H = A/2A ~ , (26) 

P = <p/<p ~ const t ^ 1 , (27) 

4 - 3r n 2 1 

3(2-r)2(2-n)hT0 • { } 

The alternative behaviour 

Q~\/(<j)-l) £ (29) 



with | < e < 2 was studied by Serna and Alimi [IT] who analyzed numerically the solutions. 

An important point to note is that it is not very meaningful to look for solutions (either 
analytically [|T(J or numerically |T7J) for Q behaving as eq. (^) or eq. ( [29] ) at all times as 



the "true" function fl may be significantly different from eq. (^5|) or eq. (|29D away from the 
basin of attraction. 



IV. THE THIRD-ORDER SOLUTIONS AND THEIR APPROACH TO GENERAL 

RELATIVITY 

The first-order being under control we now construct the third-order solutions of the 
long-wavelength scheme describing large-scale inhomogeneities. We shall set S*j = which 
in any case is negligible at late time. 
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In accordance with the method presented in |8j the third-order expansions for the 3- 
metric, the Brans-Dicke field, the energy density and the three-velocity are taken to be 



Hj = A(t) 



hij + f2(t)Rij + -(g 2 {t) - f 2 (t))Rhi 



(30) 



= 0(t) [1 + <fo(t)R] , p = pit) + p 2 (t)R , u % = u 3 (t)ViR , (31) 

where Rij is the Ricci tensor formed from the seed hy and R the corresponding Ricci scalar 
and where A(t), 0(i) and p(t) are the (Friedmann) first-order solutions of the previous 
section. 

The procedure is to insert the above expansion in eqs. (f7|)- (^0|) and keep only terms with 
two spatial derivatives. Gathering all "like" terms, i.e., the coefficients of R\ and R, gives 
the following ordinary differential equations for f 2 , g 2 , and 02 '■ 

f 2 + (3H + P)f 2 = -j , (32) 

p /2 — 3r\ 

g 2 + (2H+-)g 2 = n( — X - [30 2 + (2fiP + 3H)<p 2 + Pfi0 2 ] , (33) 



2 /a V4-3r 

kp 9 / 2 — r \ 

g 2 + (6H + —)g 2 + j = ^ (j—^ ) X - 3[0 2 + (2P + 5i?)0 2 ] , (34) 

with 

X = <j> 2 + I 2P + 3H + - J 2 + — I P + 2P 2 + PR + J 2 + -Pg 2 . (35) 

Eq. (^) can be integrated rightaway. Unfortunately the solutions to the other two equa- 
tions are not so easily forthcoming. However the interesting regime is when Q —>■ 00 (with 
(dQ / d<p) / Q 3 —>■ 0) for t big which is the limit where the first-order solutions approach their 
general relativistic limit. Now from eqs. ( |T7| - |T%| ) it follows that 3(2 — Y)H/P = 2(g f — rj)/rj 
and hence diverges as rj — ► 00 when Q — > 00, that is when g increases faster than rf. In this 
limit eqs. ([33 H341) decouple, yielding 

2 — 3r 

g 2 + 3THg 2 ~ . (36) 

As for eq. (^) it reduces to 

h + 3^/ 2 - -2/ A . (37) 



Eqs. (PB|-|37D are the same as those obtained by Comer et al. (eq. 2.15 in ||) in their study 
of inhomogeneities in general relativity. Thus, f 2 and g 2 ultimately have the same behaviour 
as the general relativity solutions discussed in |§. Unfortunately, it is difficult to make 
precise statements about 2 (and hence p 2 and it 3 ) even in the limit of diverging Q. We 
shall therefore turn to particular coupling functions. 
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A. Brans-Dicke coupling 



The late time behaviour of all Friedmann solutions in the Brans-Dicke theory of gravity 
is given by eqs. (p^-p3|). Eqs. ( p2] - |35|) for the evolution of inhomogeneities can then be easily 
integrated yielding 

= + ^ (38) 



92 



fa ft(2 - 3r) - 3(4 - 3r) 



4fi(l-/i)(3-2/i) 



2(1— ft.) 



(39) 



02 = 7l + 72^" 



3h—q 



+ 



4-3r 



02 



3[n(2 - r) - (4 - 3r)] t 



(4 - 3r)t 2 ( 1 - fc ) 

"4Q(3- 2h)(l - h) 



(4 - 3r) 2 - tt(3T 2 - 12r + 4) 

(4-3r) 2 + fi(2-r)(2 + 3r) 



(40) 



where cui, Q2, /?2, 7i, 72 are integration constants. One can check that in the limit 
Q — > oo the results for /2 and (72 agree with those of Comer et al. ||. Furthermore, when 
T = 1, corresponding to a dust fluid, all three agree with the solution of Soda et al. [[TJJ. 

The growth or decay of f'2, gi and 02 depends on the signs of the exponents in eqs. 
( ^8| - |40|) . Recalling the condition of positivity of A (eq. (0)) one finds that (/2, #2, ^2) all 
decay as time increases if and only if 



T < 2/3 , 



(41) 



that is if matter is of an inflationary type (the same condition for decay holds in general 
relativity, see |§), and if 



n > 



(4 - 3Tf 



(2-r)(2-3T) 



(42) 



This condition on Q is not demanding if p = — p (T = 0) but becomes all the more stringent 
as T -> 2/3. 

As for the three- velocity its time-dependence is obtained from 



87rrp 



1 



1V 2 



H + P 



i-r 



(43) 



which yields u 3 oc t 1 +q+H 3r ~ 2 ) — > £ 3_4 / 3r when f2 — ► 00 in accordance with the general 
relativity results ||, unless T = 1 (dust) in which case M3 = (since dust can be comoving 
in a synchronous gauge). Finally the time behaviour of the energy density correction p2 is 
given by P2 oc t q+h ^ T ~ 2 ^ — > t 2 ( 1_2 / 3r ) ; again in agreement with the general relativity result. 
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B. Dynamical coupling: an example 

Let us consider as an example the dynamical coupling generated by the function g = Erf 1 
T6|,p| . The late time behaviour of the long- wavelength, Friedmann solution, is given in eqs. 
^6| - |28|) . As already mentionned, as Q — > oo in this limit, P/H — > so that the functions fi 



and #2 are governed by the same equation as in general relativity, eqs. ( pq - pTD , the general 
solutions of which are 

or 2 

h = «x + « 2 t- (2 - r)/r - ^zr/-* /3T (44) 

and 

* = A+A«-'-j^^. (45) 
As for the equation for 2 [eqs. ( |33] - |35|) 1 it can be integrated and yields 

-4(n-2)(2-D -2(2-r) (4 - 3T) [2(2 - T) - Tl]^" 2(2 "rn"" 2) - 1 

02 = 7i* 3r " + 72^ 3r - 



6(2-r)[(2r-3)n + 2(2-T)] 



9r 2 n(4-3r)[3r(2-r) + n(3r-2)] -^-^-m-p] t 

-t 3nr . (46) 



4(9r - 4) (2 - T)(3r - 2)[(9r - 10)n + 12(2 - r)] 
All terms decay if T < 2/3, that is for inflationary matter. However 2 can still decay if 

l<r<f^f<4/3. 

V. APPENDIX: LONG- WAVELENGTH SCHEME IN THE EINSTEIN FRAME 

In the Jordan frame, the action is given by 

S= lkJ d ' x 'f- Zg " <TM0)^<R0] + S m . (47) 

In this appendix, contrarily to the main text, all the physical quantities will be labelled by 
a tilde, whereas the non-physical quantities will remain bare. Upon using a change of the 
metric and of the fields defined by 

~ 9lw = e^g^ , (48) 
(j)- 1 = e 2aM (49) 

and 

a 2 = (2w + 3)' 1 , (50) 
with ct((p) = da/ dip, one can reexpress the above action in the form 

S =^J d A x^ [R - 2gTd^dM + S M [i/> m , e 2a ^g^} . (51) 

9 



This form is said to correspond to the "Einstein frame" because the curvature term in this 
new action is the same as in the Einstein-Hilbert action for general relativity. The advantage 
of this formulation is that most of the equations can be obtained almost directly from the 
equations of |§ by combining the results for a perfect fluid and for a scalar field. 
The Brans-Dicke theory corresponds to a "linear" coupling in the sense that 



a(ip) = aip 



(52) 



Then 



-2a(p 



The field equations in the Einstein frame read 



and 



(53) 



(54) 



(55) 



with k = 8tc. In the case of a perfect fluid, this gives 



Rpv = 2d^(pd u ip + n(p + p)u^u v + ~{p - p)g (lv , 



(56) 



(57) 



One must however be aware that the p and p that appear here, defined by T^ v = (p + 



p)u^u v + pg^, with g^ u u 



•1, are not the physical energy density and pressure, but 



are related to the latter by p = e 4a p, p = e 4a p. 

In a synchronous gauge (adopting all the notations of the main text but for the unphysical 
metric) 



+ ^j=_j t (^Kf) = ^p {2Tui Ul + S]{2- r)] + 2d^d^ 



(58) 



-K + -K { K j 
2 T 4 3 1 



up 



T — TukU 



- 2<p 2 



(59) 



- (K^ - = KpTui\J 1 + u k u k - 2ipditp 



(60) 



and 



(p + l -K(p - V.VV = ^a(3p - p) 



(61) 



At first-order the trace of eq. (|58|) yields 
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d ( A\ 3r / A\ 2 2-r 



A ' V2 Ft [ A3/2 j) = Kp(2 ~ r) • (62) 

Using this equation, eq. fl59|) at first-order then gives the following equation for A: 

at \aj + T UJ + ^^" 3 + < 2 - r >^ 2 = o • (63) 

In contrast with pure general relativity, we now have an additional term coming from the 
scalar field. One thus needs an extra equation to close the system and solve for A. This is 
of course given by eq. (|6T| ) , which reads at first-order 

(p + 3Hip = ^-a(3T - 4)p . (64) 
Combining with eq. (0), one finds 



Let us now consider the Brans-Dicke case, i.e., when a is a constant. Then the previous 
equation can be integrated immediately, so as to give 

A 2(2 -T) . C . 

A-^r^f^A^- {) 

In the simplest case C = and S^Sj = 0, one finds 

A - 2Ar 1 A - 2(2 ~ r) (67) 

A~ 2Xt ' A -3r(2-r) + ^ (3 r-4)2 ' (67) 

• , i , 2a(3r-4) 

" = ft • t= 3T(2-r) + <»'(3r-4)» - (68) 



Let us make explicit the relation between this result and eqs. ( p2| - |23D in the main body 
of this article. The relation between the "Jordan time" t and the "Einstein time" t is given 
by 

i ~ t Q?+1 . (69) 

Hence 

A ~ t^e+r" , (70) 
which agrees with eq. (p2|), using a 2 = Finally, 



2aj 

~t , (71) 



which agrees with eq. 
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The third-order is obtained by taking into account the terms neglected previously, in 
particular the Ricci tensor of the three- metric. It is consistent with the approximation to 
take the Ricci tensor of in place of ^Rj, and convenient to write it in terms of the 
Ricci tensor (denoted Rij) of hy, seen as the components of a metric, to be called "seed" 
metric. We will thus look for corrections to the metric of the form 

^ = A(t) |/ a (*)i2y + i [g 2 (t) - f 2 (t)] Rh^ , (72) 

where R\ = h lk Rkj and R = R\ and corrections to the matter of the form 

<p = <p(t) + <p 2 (t)R , p = p(t) + p 2 (t)R . (73) 

We substitute 7^ = ^'7^ + ®7y in the Einstein equation and keep the terms with two 
gradients. The traceless part of eq. (ESI) gives an equation for f 2 , 



Eq. (|59"D then yields 



h + \\k - -\ • (74) 



Kp 2 = - _ 3r (g 2 + 2Hg 2 + 80O02) • (75) 



Using this relation, the trace part of eq. (EH) gives an equation involving g 2 and ip 2 : 



h + \T^2 + 6(2 - T)4m = 2 ~l^-- (76) 



Eq. (pT]) provides, once more after use of eq. ([75]), another such equation: 

(77) 



1 3r - 4 



^2 + ^^2 + 4002 



In the case of Brans-Dicke the system of equations can be solved explicitly. All the 
particular solutions for f 2 , g 2 and <p 2 behave like t 2 ~ 2X . After transformation of the time 
according to eq. (0), one sees that this is in agreement with the results in the "Jordan" 
frame. 
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